function retZ = rotZ(theta)

retZ =  [cos(theta)       -sin(theta)     0; 
         sin(theta)        cos(theta) -sin(theta); 
          0                 0            1];
end